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JL  he  current  and  future  geostationary  operational 
environmental  satellites  (GOES)  of  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA) 
are  designed  to  produce  visible  images  of  the  earth 
with  linear  arrays  of  eight  detectors.  Because  the 
imaging  instruments  are  not  calibrated  radiometri- 
cally  in  orbit,  differences  among  instrument  gains 
associated  with  the  different  detectors  may  cause 
artifcial  stripes  to  appear  in  the  images.  In  the 
data  processing  on  the  ground,  the  images  are 
“normalized  ”  to  remove  the  stripes.  Images  from 
future  geostationary  satellites,  GOES  I-M,  will  he 
normalized  by  the  method  of  matching  empirical 
distribution  functions  (EDFs).  In  this  paper  we 
report  on  a  study  of  EDF  matching  with  data  from 
GOES-7.  The  technique  was  used  to  generate  a 
normalization  look-up  table  from  data  taken  on 
^18  May  1988,  and  the  table  was  applied  to  image 
data  obtained  2  weeks  later,  on  I  June  1988.  This 
removed  the  .stripes  from  the  image.  The  technique 
is  expected  to  he  even  more  effective  with  data 
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from  GOES  I-M  because  of  improvements  in  in¬ 
strumentation. 


INTRODUCTION 

Both  the  current  and  future  Geostationary  Opera¬ 
tional  Environmental  Satellites  (GOES)  of  the  Na¬ 
tional  Oceanic  and  Atmospheric  Administration 
(NOAA)  are  designed  to  carry  instruments  that 
image  the  full  disk  of  the  earth,  or  sections  of  it,  in 
the  visible  region  of  the  electromagnetic  spectrum. 
The  current  satellite,  GOES-7  (Ensor,  1978),  spins 
at  100  rpm  on  a  north-south  axis.  Visible  radiation 
is  detected  with  an  array  of  eight  photomultipliers 
in  a  north -south  line.  The  field  of  view  from  each 
photomultiplier  is  approximately  0.8  km  sq  on  the 
earth  at  the  nadir.  With  each  rotation  of  the 
satellite,  eight  parallel  adjacent  lines  are  swept  out 
in  the  west-to-east  direction.  After  each  rotation  a 
scan  mirror  is  stepped  to  displace  the  fields  of  view 
in  the  north -south  direction.  In  the  early  1990s, 
the  GOES  I-M  system  (Komajda  and  McKenzie, 
1987)  is  expected  to  become  operational.  Instead 
of  spinning,  those  satellites  will  be  three-axis  stabi¬ 
lized  and  will  always  present  the  same  side  to  the 
earth.  Each  GOES  I-M  imager  will  utilize  a  linear 
array  of  eight  detectors,  which  will  be  silicon 
photodiodes  instead  of  photomidtipliers,  and  it  will 
scan  both  west-to-east  and  east-to-west.  Despite 


90 


0(m-42^7/H9 /$<).<)() 

This  fjrtiric  bs  in  fhc  Pxddir  Dotudin. 


185 


186  Weinreh  et  al. 


MkmxOJhni  (GOES  a-?) 

IJihniililhni  (OOeS  l-M) 

Figure  1.  Scanning  geometry  of  GOES  imagers  (not  to 
scale). 

such  differences  in  instniinentation,  the  current 
and  future  GOES  imagers  will  have  similar  scan¬ 
ning  geometries,  depicted  in  Fig.  1. 

Henceforth,  we  will  refer  to  each  detector  and 
its  associated  optics  and  signal-processing  electron¬ 
ics  as  a  channel.  Each  imaging  instrument  is  then 
said  to  have  eight  visible  channels.  When  stimu¬ 
lated  with  incident  radiation,  each  channel  re¬ 
sponds  with  an  output,  measured  in  digital  counts. 
W^e  use  the  term  “gain”  to  refer  to  the  ratio  of  the 
change  in  the  coimt  output  to  the  change  in 
intensity  of  the  incident  radiation.  The  visible 
channels  of  the  GOES  I-M  imagers  are  designed  to 
be  linear,  i.e.,  the  output  will  be  a  linear  function 
of  incident  intensity,  and  the  gain  will  be  indepen¬ 
dent  of  intensity.  However,  the  visible  channels  of 
GOES-7  were  designed  with  the  output  approxi¬ 
mately  equal  to  the  square  root  of  the  incident 
intensity  to  take  advantage  of  the  increased  sensi¬ 
tivity  of  the  photomultipliers  at  low  intensity  levels 
(Ensor,  1978;  Bristor,  1975).  The  photomultipliers, 
which  typically  operate  in  a  slightly  nonlinear  fash¬ 
ion,  may  also  add  to  the  nonlinearity. 

Once  a  GOES  is  in  orbit,  its  visible  channels 
cannot  be  calibrated  radiometrically,  because  there 
is  no  calibrated  source  of  visible  radiation  on  board. 
Since  the  gains  (and  offsets)  of  the  eight  channels 
may  not  all  be  equal  and  may  change  with  time, 
artificial  east -west  stripes  can  appear  in  the  im¬ 
ages.  Figure  2  shows  striping  in  an  image  from 


GOES-7.  Striping  may  also  affect  image  data  from 
GOES  I-M,  although  it  is  expected  to  be  less 
severe.  One  reason  for  this  is  that  the  gains  of  the 
photodiode  detectors  of  GOES  I-M  are  expected  to 
be  more  nearly  uniform  and  more  stable  in  time 
than  the  gains  of  the  photomultipliers  of  CiOES-7. 
Also,  the  data  from  CiOES  I-M  will  be  digitized 
more  finely.  The  6-bit  word  stnicture  of  GOES-7 
will  be  increased  to  10  bits  on  GOES  I-M.  The 
more  bits  per  word,  the  less  significant  the  striping 
from  quantization  error. 

In  operational  data  processing  at  NO.\A's  Na¬ 
tional  Environmental  Satellite,  Data,  and  Infonna- 
tion  Service  (NESDIS),  we  compensate  for  chan- 
nel-to-channel  differences  in  gain  with  a  procedure 
called  nonnalization,  which  is  applied  to  the  data 
with  look-up  tables.  It  is  done  automatically,  con¬ 
tinuously,  and  in  real  time  at  the  ground  station  at 
Wallops,  VA,  as  the  data  are  received  from  the 
satellites  and  then  retransmitted  to  users.  Lx)ok-up 
tables  are  generated  off  line  (and  not  in  real  time) 
in  the  main-frame  computers  at  Suitland,  MD. 
For  each  .satellite  instalment,  a  single  table  is 
applied  to  all  image  data  from  the  entire  disk  of 
the  earth.  Because  characteristics  of  the  instni- 
ments  change  with  time,  new  look-up  tables  are 
generated  occasionally,  usually  once  a  week  or  less 
often.  The  current  method  of  generating  the  oper¬ 
ational  look-up  tables  (J.  Lienesch,  NOAA/ 
NESDIS,  personal  communication,  1980)  has  been 
in  use  since  the  mid  1970s.  However,  NESDIS 
plans  to  apply  a  different  method,  namely,  the 
matching  of  empirical  distribution  functions 
(EDFs),  to  generate  the  tables  for  GOES  I-M, 
because  it  is  more  efficient. 

Matching  of  EDFs  is  a  standard  statistical 
technique  and  has  already  been  applied  to  satellite 
data  from  the  Landsat  Multi.spectral  Scanner  (MSS) 
[see,  e.g.,  Horn  and  Woodham  (1979)].  In  that 
work,  a  normalization  look-up  table  is  derived  from 
and  applied  to  the  same  sector  of  image  data.  Our 
application  is  different,  however,  because  we  are 
required  to  derive  a  table  from  a  single  “depen¬ 
dent"  sector  on  a  particular  day  and  apply  it  for 
weeks  afterw'ards  to  independent  image  data  from 
all  over  the  earth's  disk. 

In  preparation  for  the  launch  of  CiOES  I.  we 
have  been  studying  the  feasibility  of  the  EDF 
technique  by  applying  it  to  data  from  GOES-7, 
which  was  the  only  multichannel  satellite  imager 
who.se  data  were  conveniently  available  to  us.  Data 


Figure  2.  Unnormalized  GOES-7  image.  18  May  1988. 


from  GOES-7  offer  a  severe  test  of  the  method,  assumption,  the  basic  premise  becomes  that  the 
because  it  is  expected  that  the  GOES  I-M  sensors  distributions  of  the  outputs  of  each  channel  should 
will  be  better  behaved  and  more  stable  in  time  be  identical. 

than  those  of  GOES-7,  as  we  explained  previously.  In  our  approach,  we  designate  one  of  the 

This  paper  summarizes  some  of  the  results  from  channels  as  a  reference  channel.  Then  the  outputs 

the  GOES-7  study.  More  details  can  be  found  in  of  the  other  channels  are  adjusted  with  the  normal- 

Weinreb  et  al.  (1989).  ization  tables  so  that  their  distiibutions  are  the 

same  as  that  of  the  reference  channel.  The  refer¬ 
ence  channel  should  be  selected  on  the  basis  of  its 
relative  stability,  low  noise,  and  maximal  use  of  the 
THEORY  dynamic  range  of  the  data  system  without  clipping 

at  either  the  low'  or  high  ends. 

The  basic  premise  is  that  if  several  channels  view  To  generate  a  nonnali/ation  l(M)k-up  table,  we 

the  same  scene,  their  outputs  should  be  equal,  and  |)egin  by  selecting  a  sample  of  full-resolution  un- 
this  should  be  the  ca.se  regardless  of  how  bright  the  normali/etl  earth-scene  data  covering  as  much  of 

scene  is.  In  actual  application,  no  two  channels  the  range  of  intensities  as  possible.  For  GOES  I-M. 

ever  view  the  same  scene.  Instead,  we  assume  that  the  area  will  be  rectangular,  extenditig  several 

with  a  large  ensemble  of  measurements,  the  distri-  thousand  pixels  both  east  to  west  and  north  to 

bution  of  the  intensity  of  the  earth  radiation  inci-  south.  Corresponding  to  the  incoming  radiance 

dent  on  each  detector  will  be  similar  (Horn  and  frotn  any  pixel,  the  instniment  will  resijond  with 

Woodham,  1979).  (In  practice,  the  distributions  an  output  .x,  in  digital  counts.  One  can  compile 

will  not  be  identical,  but  the  larger  the  ensemble,  the  discrete  density  function,  i.e..  the  histogram, 

the  more  similar  they  will  become.)  With  that  describing  the  relative  fre(inency  of  oct  iirrence  of 
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each  possible  count  value,  for  each  channel.  For 
Channel  i,  which  is  the  channel  to  be  normalized, 
let  the  histogram  be  p,(x).  An  empirical  distribu¬ 
tion  function  (EDF)  P,(x)  can  then  be  generated; 
viz., 

^’,(x)=  t  p,(0. 

(  =  0 

(The  EDF  is  also  known  as  a  cumulative  histogram 
of  relative  frequency.)  The  EDF  is  a  nondecreas- 
ing  function  of  x,  and  its  maximum  value  is  unity. 
For  convenience,  however,  we  have  chosen  the 
maximum  value  to  be  100%;  i.e.,  if  the  maximum 
possible  output  in  counts  is  X,  then 

P,(X)  =  100%. 

In  these  terms,  the  basic  premise  says  that  for 
each  output  value  x  in  Channel  i,  the  normalized 
value  x'  should  satisfy 

P.(.t')  =  P,(x),  (1) 

where  the  subscript  r  refers  to  the  reference  chan¬ 
nel.  In  practice,  not  only  is  P,  nondecreasing,  but 
it  is  also  monotonically  increasing  as  a  function  of 
x'  in  the  domain  of  x'  where  there  are  data. 
Therefore,  it  can  be  inverted,  yielding  the  solution 
for  x'. 

x'=P,  ‘(P,(x)).  (2) 

WTien  it  is  applied  sequentially  for  every  possi¬ 
ble  count  value  x,  Eq.  (2)  generates  the  normaliza¬ 
tion  look-up  table  relating  each  x  to  an  x'.  Figure 
•3  depicts  how  the  procedure  is  applied  in  actual 
practice  to  generate  one  entry  in  the  table.  The 

f  igure  3.  Illustration  of  procedure  to  generate  nonualization 
look-up  table  (see  text  for  explanation) 


P  (%) 


figure  depicts  idealized  EDFs  for  the  reference 
channel  and  Channel  i.  In  the  figure  the  EDFs  are 
continuous,  but  in  practice  they  are  discrete,  being 
specified  only  at  integer  values  of  x.  To  find  x,'. 
the  nonnalized  count  value  corresponding  to  the 
observed  count  value  of  Xj,  the  following  is  the 
procedure:  First,  for  the  count  value  X[  in  Chan¬ 
nel  i,  find  the  percentage  vahie  from  the  EDF  of 
Channel  t.  In  the  illustration  it  is  P,(  Xi)-  Then  find 
the  point  on  the  reference  channel’s  EDF  with  the 
same  percentage  value.  According  to  Eq.  ( 1 ),  that 
percentage  can  also  be  expressed  as  Prlx/).  Fi¬ 
nally,  use  the  EDF  of  the  reference  channel  to  find 
the  normalized  count  value  x,'.  Since  the  data  are 
actually  discrete,  we  will  need  to  interpolate  within 
the  EDF  of  the  reference  channel  to  find  the  value 
of  Xj',  which  must  then  be  rounded  to  the  nearest 
integer. 


RESULTS  AND  DISCUSSION 

A  normalization  look-up  table  was  generated  from 
a  dependent  sample  consisting  of  unnormalized 
data  from  GOES-7  on  18  May  1988.  The  table  was 


f  igure  I  Location  on  glolse  of  dependent  and  independent 
samples  tif  image  data.  1996  X  24f)0  pixels. 
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Generation  of  Normalization  I^k-Up  Table 


Figure  6.  Empirical  distribution  functions  for  unnonnalized 
GOES-7  image  data,  18  May  1988. 


applied  to  normalize  an  independent  sample  of 
image  data  from  GOES-7  2  weeks  later,  on  1  June 
1988.  Both  samples  were  produced  at  the  same 
time  of  day,  approximately  1400  EDT  (1800Z), 
when  the  disk  of  the  earth  beneath  the  satellite 
was  in  full  sunlight. 


The  dependent  sample  consists  of  unnormalized 
data  from  all  pixels  in  a  rectangular  sector  1996 
pixels  east -west  and  2400  pixels  north  south,  de¬ 
picted  at  the  equator  in  Fig.  4.  An  image  of 
approximately  one  sixth  of  that  sector  appears  in 
Fig.  2,  and,  as  we  saw  previously,  it  is  striped.  The 
striping  is  most  severe  at  the  mid-  and  high-inten¬ 
sity  levels  (low  and  high  cloud),  which  appear  gray 
and  white,  respectively. 

For  each  of  the  channels,  histograms  were 
compiled  from  the  data  in  the  full  dependent 
sample.  Each  histogram  represents  data  from 
598,800  pixels.  The  three  panels  in  Fig.  5  show 
the  histograms  for  Channels  2,  5,  and  6.  (The 
abscissae,  labelled  “Intensity  (Counts),"  are  the 
output  levels.)  Based  on  the  criteria  mentioned 
previously,  Channel  2  was  designated  as  the  refer¬ 
ence.  We  chose  Channels  .5  and  6  as  examples 
because  their  histograms  are  the  most  different 
from  the  reference  channel’s.  The  histogram  of 
Channel  5  is  displaced  towards  the  high  intensities 


Figure  7.  Unnonnalized  GOES-7  image,  1  June  1988. 
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relative  to  that  of  the  reference  channel,  its  upper 
end  is  clipped,  and  it  is  broader  than  the  reference 
channel  s,  indicating  that  the  gain  in  Channel  5  is 
greater  overall  than  that  of  the  reference  channel. 
With  rhannel  6,  the  situation  is  ••eversed,  because 
the  gain  is  lower  overall  than  that  in  the  reference 
channel.  The  diffeiC'ices  among  the  three  his¬ 
tograms  also  suggest  that  the  dependence  of  gain 
on  intensity  (or  the  degree  of  nonlinearity)  varies 
from  channel  to  channel. 

Figure  6  shows  the  EDFs  of  the  three  chan¬ 
nels,  which  were  computed  from  the  histograms  in 
Fig.  5.  The  abscissa  is  the  output  level,  and  the 
ordinate  is  the  percentage  of  the  data  with  outputs 
at  or  below  that  level.  We  applied  Eq.  (2)  to  these 
EDFs  and  those  of  the  other  channels  to  generate 
the  nonnalization  look-up  table.  [See  Weinreb 
et  al.  ( 1989)  for  a  listing  of  the  table.]  That  table 
will  be  applied  to  independent  image  data,  as 
described  below. 


Application  to  Independent  Sample 

The  independent  sample  of  unnormalized  GOES-7 
data  was  produced  on  1  June  1988,  from  the  upper 
sector  of  1996  x  2490  pixels  depicted  in  Fig.  4.  An 
image  of  approximately  one  sixth  of  that  .sector  is 
shown  in  Fig.  7.  As  expected,  there  is  striping, 
and  it  is  most  severe  at  the  mid-  and  high-intensity 
levels.  Figvires  8  and  9  show  the  histograms  and 
EDFs,  respectively,  compiled  in  Channels  2,  5, 
and  6.  Their  shapes  are  quite  unlike  the  corre- 


Figure  9.  Empirical  distribution  functions  for  unnonualized 
GOE.S-7  image  data,  1  June  1988. 


sponding  histograms  and  EDFs  from  the  depen¬ 
dent  sector,  reflecting  the  differences  in  the  cloud 
distributions  in  the  two  sectors.  However,  for  the 
nonnalization  look-up  table  to  be  effective  on  the 
independent  sector,  the  chamiel-tOK'hannel  rela¬ 
tionships  among  the  EDFs  must  be  similar  for  the 
two  .sectors;  i.e.,  the  relative  (channel  to  channel) 
gain  functions  must  be  similar,  .\lthough  it  may 
not  be  obvious  at  this  point,  the  results  below 
demonstrate  that  this  is  in  fact  the  case. 

We  nonnalized  the  1  June  data  by  applying 
the  18  May  look-up  table.  Figure  10  shows  the 
histograms  of  the  normalized  data  for  Channels  .5 
and  6.  In  position  and  shape  they  are  now  much 
more  like  the  histogram  of  the  reference  channel 
(reproduced  in  the  upper  panel),  as  they  should  be 
if  the  nonnalization  is  to  be  successful.  .An  unusual 
feature  in  these  histograms  is  the  presence  of 
breaks,  i.e.,  intensity  levels  with  a  zero  frequency 
of  occurrence.  Because  the  data  are  discrete,  breaks 
will  (Kcur  when  normalization  expands  a  region  of 
a  histogran),  as  is  explained  more  fully  in  Weinreb 
et  al.  ( 1989). 

Figure  11  shows  the  EDFs  of  the  nonnalized 
data  for  all  three  channels.  .As  it  should,  the  nor¬ 
malization  process  made  the  differences  among 
them  practically  insignificant.  (The  EDFs  of  the 
other  five  channels  behave  similarly.)  The  largest 
differences  occur  where  the  EDFs  of  Channels  5 
and  6  have  flat  spots.  These  are  caused  by  the 
breaks  in  the  histograms  and  are  an  artifact  of 
digitization. 

Figure  12  is  the  nonnalized  image  of  the  same 
area  as  was  shown  in  Fig.  7.  It  is  the  “after  to  the 
“before  of  Fig.  7.  The  improvement  is  obvious, 
since  w'e  cannot  see  any  stripes  in  the  image. 

The  effectiveness  of  the  nonnalization  means 
that  the  channel-to-channel  relationships  among 
the  EDFs  remained  essentially  the  same  between 
18  May  and  1  June,  as  was  surmised  earlier,  even 
though  the  cloud  patterns  and  the  EDFs  them- 
.selves  changed.  Therefore,  the  relative  gains  and 
offsets  amt)ng  the  channels  also  must  have  re¬ 
mained  the  same.  Ftirthermore,  although  we  do 
not  show  the  results  here,  we  found  that  the 
effectiveness  of  the  18  May  normalization  look-up 
tables  decreased  only  slightly  with  time  over  a 
6-week  period.  Therefore,  any  changes  in  relative 
gains  and  offsets,  which,  e  g.,  might  have  resulted 
from  seasonal  variation  of  temperatures  on  the 
satellite,  had  to  be  small  over  that  pericKl. 
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Figure  10.  Histograms  of  nonnalized 
GOES-7  image  data,  1  Jime  1988:  (-)  Ghan- 
nel  2;  (b)  Channel  5;  (c)  ('hannel  6.  (Data  in 
Channel  2,  the  reference  channel,  are  neser 
nonnalized.) 
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Figure  11.  Empirical  distribution  functions  for  nonnalized 
GOES-7  data,  1  June  1988:  (-)  Channel  2  (reference);  (  -  -  - ) 
Channel  5;  (---)  Channel  6. 


CONCLUSION 

The  case  study  presented  in  this  paper  is  strong 
evidence  that  normalization  by  EDF  matching  is 
an  effective  method  for  removing  striping  from 
visible  images  from  GOES.  Application  of  a  nor¬ 
malization  look-up  table  generated  from  GOES-7 
data  of  18  May  1988  removed  the  stripes  from  an 
image  obtained  on  1  June  1988.  The  method 
worked  despite  the  nonlinearities  in  the  outputs  of 
the  GOES-7  visible  channels  and  the  substantial 
channel-to<'hannel  nonunifonnities  in  gain.  We  ex¬ 
pect  it  to  work  at  least  as  well  with  data  from 
GOES  1-M,  becau.se  the  responses  in  the  GOES 
I-M  channels  are  expected  to  be  linear  and  more 
nearly  uniform  than  those  of  GOES-7.  Further,  we 
found  that  the  same  normalization  table  remained 
effective  for  a  least  6  weeks.  Since  the  gains  of  the 
GOES  I-M  channels  are  expected  to  be  more 
stable  in  time  than  those  of  GOES-7,  we  would 


Figure  12.  Nonnalized  GOES-7  image,  1  June  1988. 
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expect  a  norniali/atioii  tal)le  to  remain  valid  for 
even  loii^or  periods  uith  CX)ES  I-M. 

The  (titlhors  ure  iuclchted  to  Dung  llan  of  \OAA  /  .V/'-SD/S /or 
pnn  iding  the  rate  image  data-  to  C.ene  Leejg  of  SOA.\  / 
\ESDIS  for  producing  the  " sectorizer"  images  that  appear  in 
Figs  r.  and  12:  to  Donald  Mack  of  Integra!  Systtans,  Inc.. 
of  Lanham,  MD,  for  supplying  information  on  the  application 
of  F.DF  matching  to  data  from  Landsat;  and  to  George  Jones 
and  Paige  Bridges  of  XOAA/ \TSDIS  for  preparing  the 
figures. 
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